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Abstract 

Polaron binding energy and effective mass are calculated for semiconductors 
with wurtzite crystalline structure from the first order electron-phonon corrections 
to the self-energy. A recently introduced Frohlich-like electron-phonon interaction 
Hamiltonian which accounts for the LO and TO polarizations mixing due to the 
anisotropy is used in the calculation. The polaronic damping rates are evaluated 
for finite temperature. Numerical results are reported for GaN . It is shown that the 
electron-phonon coupling is strong enough to justify the necessity of the inclusion 
of second-order corrections. 
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Electronics and optoelectronics based on III-V wide band gap nitride semiconducting 
systems have emerged as a matter of particular interest in recent years. An increasing 
number of papers devoted to the study of these materials and their applications can be 
found in the literature (see, for example, Refs. [1 — 4]). 

Nitrides have the wurtzite as their natural crystalline structure, and exhibit highly unusual 
properties among the III-V group of semiconductors, more resembling, in some aspects, 
the II- VI compounds. It is possible to mention, for instance, the large magnitudes of the 
spontaneous and piezoelectric polarizations [5]. 

An aspect which is only beginning to be investigated in crystals with such a structure is 
the carrier-polar-optical-phonon interaction mechanism. A Frohlich-type electron-phonon 
Hamiltonian was recently put forward [6] within the non-retarded Born-Huang dielectric 
continuum approach. In that work only the three optical-phonon branches which are 
infrared- active in wurtzites are considered. Probably, the most important result arising 
from that work is the mixing of the longitudinal and transverse optical modes due to the 
crystalline anisotropy. As a consequence, both kind of modes do contribute to the carrier 
scattering, although LO-like contributions should be dominant. A similar behaviour is 
well known in semiconducting heterostructures, where the presence of interfaces breaks 
the isotropy in one or more spatial directions, leading to the mixing of the longitudinal 
and transverse oscillation polarizations [7, 8] , and all of them have to be included in the 
study of the electron-phonon interaction because even the TO-type modes can provide a 
non-negligible contribution [9] to the scattering rates and polaron properties. 

Furthermore, the existence of a spatial direction with lesser symmetry also gives rise to 
phonon dispersion. To illustrate, let us mention the application of the macroscopic dielec- 
tric continuum model to the study of the long-wavelength optical oscillations in dielectric 
or semiconducting heterolayers [10]. In that case, even when no phonon dispersion is 
considered for the bulk constituent materials, a dispersion law w(q) appears for the so- 
called interface modes. For wurzite III-V nitrides, phonon dispersion effects should be 
expected. As already discussed in Ref. [6], the dependence of the oscillation frequencies 
upon the phonon wavevector q is not a function of the wavevector magnitude but of the 
angle between q and the c axis of the structure (referred to from now on as z direction), 
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the direction of less symmetry in the unit cell. 

The electron-phonon scattering rates [6] are sufficiently large to recognize a rather strong 
interaction, and to motivate the study of the polaronic corrections to the energy and the 
effective mass of slow electrons in the conduction band of wurtzite nitride materials. 

The polaron is a quasiparticle state consisting of an electron and its surrounding phonon 
cloud. When electrons move through a polar material, they polarize the surrounding 
medium and couple to the self-induced polarization field [11]. This coupling slows down 
the electron as if it acquired an additional mass, and -at the same time- behaves as a 
potential well, decreasing its zero momentum energy, thus tending to localize the electron 
(which really happens when the coupling is strong enough). Even for weakly polar mate- 
rials, these corrections may be quite important and will necessarily have to be taken into 
account, for instance, in the study of optical and transport properties. The aim of this 
work is -then- to calculate the polaronic corrections in wurtzite nitride semiconducting 
materials according to the dielectric long-wavelength continuum model. In particular, 
numerical results will be presented for GaN . 

The electron-phonon interaction Hamiltonian corresponding to the infrared-active phonon 
modes of wurtzite structures can be written in general form as the sum of two contributions 
[6]: one of them corresponds to electron-LO-like-phonon interaction and the second one 
stands for the electron-transverse-optical(TO)-like-phonon interaction. The longitudinal 
and transverse electron-phonon vertices are given by[6] 



Here, l/e* ± = 1/eJ* - l/e° ± and 1/e* = 1/ef - l/e°; ef(ef) being the high-frequency di- 
electric constant perpendicular to (along) the z axis, and 6° = ^ijj\ L /ijj\] e °z = ^^Il/^I 
are the static dielectric constants. uj z l and uj±l are the LO-phonon frequencies along and 
perpendicular to the z axis respectively, and u z and ujj_ are the corresponding lattice 
dispersion frequencies. For simplicity, ef and are assumed to be identical [6]. 




2ne 2 h \ sin 2 9 cos 2 9 




(1) 
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The characteristic frequencies Vt L and VL T are also functions of 6, the angle between the 
wavevector q and the z direction. They are given by [6] 




■ = uo 2 zL cos 2 9 + uj\ l sin 2 9, 

n 2 2 • 2 a i 2 2/i 

iZ T = u; 2 sin 6* + cj_|_ cos 0, 



(2) 



Unlike isotropic materials, the TO-like vertex 



is, in general, different from zero due 



to the polarization mixing. A pure transversal mode is obtained for 9 = |. 

The electronic states are described, as usual, within the effective mass approximation: 
The electron wavevector k is the quantum number which labels the state with energy 
E = h 2 k 2 /2m* in the conduction band, which is assumed to be spherically symmetric 
with an effective mass m*. 

In order to illustrate the results of our calculations we choose GaN with frequencies and 
dielectric constants: u zL =735 cm" 1 , lu ±l =743 cm" 1 [12], =5.29, e° z =10.01 and 
e° ± =9.28. The conduction band effective mass is taken to be m* =0.20m [13], m being 
the electron bare mass. 

To deal with the entire range of 9, a Frohlich-like coupling constant ol is introduced 
through a\ = m* e 4 / (2e* ± e* z h 3 uj L ) . uo L would be an "effective longitudinal frequency". In 
our case, we choose the value ojl = &l(9 = 7r/4) to treat contributions coming from u±l 
and uo z l on equal footing. For GaN, Tiul = 91.64 meV and «l = 0.39. An "effective 
longitudinal polaron radius" pl introduced through p\ = Ti/2m*u L has the value 14.4 
lA for GaN. 

To describe the electron-TO-like-phonon interaction in terms of certain dimensionless 
coupling parameter, let us introduce another Frohlich-like constant ar, so that a\ = 
m*e 4 /(2e* ± e* z h 3 ujT)- Under the same argument used to introduce the frequency col, an 
"effective transversal frequency" uj t is chosen to be ojt = Cl T (9 = 7r/4). 

For GaN we get Tiuj t = 67.86 meV, leading to the values cct = 0.45, and p T = 16.75A, 
for the coupling constant and the "effective transversal polaron radius", respectively. 
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Now, the longitudinal and transverse electron-phonon vertices given in Eqn. (1), modify 
in the following way: the factor 2%e 2 h is changed to 4iih 2 aL((l/ 'e*)(l/ el)) -1 / 2 Pl^l and to 
47r^ 2 o;T((l/e*)(l/el))~ 1 / 2 PT^T for the longitudinal and the transversal case, respectively. 

This way of choosing uj^, out, and «y, to represent the electron-phonon interaction for 
the wurtzite is, of course, not unique. It is possible to follow different criteria; but what 
remains clear is that the final numerical results will be the same. Indeed, it makes no sense 
to compare the strengths of the electron-LO-like-phonon interaction and of the electron- 
TO-like-phonon interaction through the values of ol and ax- However, both ai and «t 
can be representative values if comparisons between the strength of the corresponding 
interactions are made between different wurtzite materials. 

Calculating the first order correction to the electron self-energy, the real and imaginary 
parts will be given by [14] 



Re 
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; (3) 



with x = cos 9, cos v = cos(</? — 7) sin 77(1 — x 2 ) 1 ^ 2 + xcosrj, is the cosinus of the angle 
between the wavevectors k and q, 7 is the polar angle associated to the electron wavevector 
k, and r\ is the angle between k and the z direction. The direction dependent characteristic 
phonon frequencies are given by 



n 2 (x) = i 



' (u 2 zL - u 2 ±L )x 2 + u 2 ±L , S = L; 



{ t^ ± -u 2 )x 2 + ^ z , 



S = T, 



(4) 
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and 



fs(x) = 



l-x 2 



*) 2 x 2 (l-x 2 ) 



I (e'i-e-V 2 x2 + ( £ 0_ e oo )a; 2 (1 _ a .2 ) , 



(5) 



The polaron binding energy and effective mass are readily calculated from the real part of 
the electron-phonon retarded self-energy in the limits of one particle and zero temperature 
[12]. As usual, only the terms up to second order in k are kept. Accordingly, the polaron 
binding energy will be given by, 



e p = -(a L £ L huj L + a T ^Thu T ), 



(6) 



On the other hand, the polaron mass is given by 



m„ = 



p 1 - Hp(v) ' 



(7) 



H P (r]) = (a L n L1 + a T /iri) cos 2 rj + (a L fi L 2 + «t Vti) sin 2 77. 



(8) 



Hence, the polaron effective mass in wurtzite crystals, with an electron-polar-optical- 
phonon interaction described by Eqn. (1) is a function of the angle between the electron 
wavevector k and the c axis of the crystal lattice. The expressions for the quantities £s 
and jis are the following (S = L,T): 



1/2 ,1 



<5£iy" r fs(x) 
u s J Jo n 3 J 2 (x) 



dx; 



(9) 



1 / * M l/2 f 1 X 2 fs{x) 



HS2 = ^(u S e* z el) 1/2 j{] 



i(l-x 2 )f s (x) 



dx; 



(10) 
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For GaN, the following numerical results were obtained: = 9.86 x 10 1 ; £t — 3.79 x 
10- 3 ; i_i L1 = 1.68 x lO^ 1 ; fi L 2 = 1.61 x 10 _1 ; Uti = 8.08 x 10~ 4 , and fi T2 = 5.46 x 
10~ 4 . Then, the numerical value of the polaron binding energy will be e p = —(35.25 + 
0.12) meV = —35.35 meV. Fig. 1 shows the relative magnitude of the effective mass first 
correction, Mr = m p /m* — 1, versus rj for GaN. It is seen that the average magnitude of 
the increasing of the electron effective mass is around 7%, and the angular variation of 
this quantity is of approximately 4%, from 9 = to 9 = n/2. 

The amount of the total corrections due to the electron-phonon interaction is, in fact, 
rather significative. Indeed, the relative weight of the contributions coming from the TO- 
like modes is very small: 0.34% to the binding energy correction and 0.38% -in average- 
to the effective mass correction. As can be expected, the main contribution comes from 
the LO-like modes. 

Anisotropy in the polaron effective mass is shown to be, in fact, a very small effect. 
It is worth noting that in this work an isotropic spherical effective mass is assumed. 
However, in wurtzite structures the electron effective mass should have different values 
for m* and rn]_. This fact could be easily included in this formalism taking into account 
in Eq.(3) the corresponding electron wave function and energy dispersion. Nevertheless, 
for GaN both masses are almost the same [15]. Therefore, the result expressed by Eq. (7) 
is valid at this level of approximation. If only the LO-phonons were active, the polaron 
effective mass correction is Mr(7r/2) ~ Mr(0) = 7.52 x 10~ 2 . By contrary for TO-phonons 
Mr(vr/2) = 3.83 x 10" 5 and Mr(0) = 5.61 x 10" 5 , showing that the interaction with TO- 
like phonons is responsible for the introduction of a strong ^-dependence. However, the 
relative contribution of this last correction makes the total result for the anisotropy to be 
small. 

The magnitude of the polaron mass anisotropy is also related to some other material 
factors like the phonon frequencies and dielectric constant. Other III-V wurtzite semicon- 
ductors like AIN and InN with similar ionicities show roughly the same results[16]. By 
contrast wurtzite II- VI semiconductors which have larger ionicities but similar phonon 
frequencies, should present a larger polaron mass anisotropy, which could be calculated 
in the same way, but experimental Frohlich constants for those materials are scarce. It 
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can be expected that for the case of different wurtzite materials, a higher influence of the 
electron-phonon interaction in the effective mass anisotropy could be appreciable. 

In the results of Ref. [6] for the scattering rates at T = 300 K, it is also seen that the 
contribution coming from the TO-like mode is comparatively small. The largest value 
of the TO-like matrix element is reported to be 7% of the LO-like matrix element. This 
indicates that for finite temperatures, the relative amount of the TO-like contribution 
might increase. To investigate that point, we have computed for k = the so-called po- 
laronic damping rate in wurtzite GaN from the first order imaginary part of the retarded 
self-energy of Eqn. (3). In our case -as can be seen in the corresponding expression-, only 
the phonon absorption process is non-vanishing because for k the emission threshold is 
never reached even for finite temperature. The polaron damping rate is given through the 
equation: 



with 13 = l/k B T. 

Fig. 2 shows the dimensionless quantity Is as a function of the temperature. Fig. 2(a) 
corresponds to the LO-like case, and Fig. 2(b) corresponds to the TO-like case. It is 
inmediately seen that the value of the TO-like contribution to the polaron damping rate 
is always much smaller than the LO-like one; but its magnitude increases in one order of 
magnitude when going from T = f 00 K to T = 300 K, and this could be an important 
effect for a different wurtzite material, where the combined values of the parameters give 
rise to a larger electron-TO-like-phonon matrix element. 

In summary, we have reported in this work the polaron binding energy and effective mass 
in first order perturbation theory for wurtzite materials. Numerical results for the case of 
GaN show that the magnitude of the coupling between the conduction electrons and the 
polar optical phonons is large indeed. Both, the LO-like and TO-like oscillation modes 
contribute, even when it is seen that the relative weight of the TO-like contribution 
is much smaller than the one coming from the LO-like phonons; but it is possible to 



I 



as^si^l) Is] 




(11) 
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expect appreciable contributions from the TO-like phonon scattering processes at room 
temperature. The value of the polaron binding energy indicates the great importance that 
the carrier-phonon interaction might have, for instance, in the study of exciton-related 
optical processes in wurtzite semiconducting materials and heterostructures based on 
them. It is very well known that these materials are the subject of interesting applications 
in optoelectronics, where III-V nitride-based blue-green semiconductor lasers have already 
been produced. 

Obviously, the electron-phonon coupling is strong enough to allow for non-negligible con- 
tributions from the second order diagrams in the self-energy expansion. According to the 
perturbation Hamiltonian [6], these second order diagrams, with two phonon lines, occur 
with three different structures: two LO-like phonon lines, two TO-like phonon lines, and 
one LO-like plus one TO-like phonon lines. For zero temperature, in the evaluation of the 
second order corrections, it is necessary to consider only the contributions coming from 
the first of these structures: that of two LO-like phonon scattering processes. This, in the 
case of GaN, is justified by the numerical value above reported for the polaron binding 
energy in first order approximation (not to mention that the magnitude of the effective 
mass correction is considerably smaller). Nevertheless, for finite temperature, the contri- 
butions of TO-like phonons in second order might also become significative, possibly for 
a different wurtzite material. The work of calculating these second order corrections, and 
their inclusion in the study of exciton-related optical absorption is already in progress. 
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Figure Captions 

Fig.l: Relative polaron effective mass for GaN as a function of 77, the angle between the 
electron wavevector and the c axis of the wurzite crystalline structure. 

Fig. 2: The dimensionless quantity Is as a function of the temperature for GaN: (a) 
LO-like phonon absorption process, and (b) TO-like phonon absorption process. 
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